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We present approximate analytical solutions to the Hamiltonian and momentum 
constraint equations, corresponding to systems composed of two black holes with 
arbitrary linear and angular momentum. The analytical nature of these initial 
data solutions makes them easier to implement in numerical evolutions than the 
traditional numerical approach of solving the elliptic equations derived from the 
Einstein constraints. 

Our work is based on descriptions of black holes in ingoing Eddington- 
Finkelstein (iEF) coordinates. This choice is motivated by the fact that, in these 
coordinates, surfaces of constant time "penetrate" the event horizon. Foliations 
that penetrate the horizon facilitate the excision of the singularity from the compu- 
tational domain. The essence of black hole excision is the removal of the singularity 
while preserving the integrity of the spacetime accessible to observers outside the 
black hole. n 

Matzner, Huq, and Shoemaker El presented a simple way to construct approxi- 
mate initial data for multiple black-hole systems, which consists essentially in super- 
posing the spatial metric corresponding to single black holes with arbitrary spins, 
positiptns, and velocities. In order to further reduce the residual errors, Marronetti 
et al. 13 proposed variation of the superposition method that preserves the simplicity 
of being analytical. Essentially, the method consists of multiplying "attenuation" 
functions into the recipe of the previous section. The purpose of the attenuation 
functions is to cancel the effects of a given hole on the neighborhood of the other 
hole, to eliminate the singularities present in the constraint equations. 

We show here that even for interestingly close separation scenarios, the atten- 
uation method errors are small, and in fact both the h and ^oo norms are smaller 
than those of the truncation error (the discretization error in the calculation) in 
simulations at currently accessible computational resolutions. The results shown 
in this paper correspond to a head-on collision (which is simpler to display), but 
very similar results are found at comparable separations for non head-on data. The 
parameters of the system are Mi = M 2 = M = 1, |aj| = |a 2 | = a = 0.5M along 
the z axis, and the holes are boosted against each other in the x direction with 
velocity v = 0.5. In order to evaluate the accuracy of these approximations, we 
calculate the residuals of the Hamiltonian and momenturmxonstraint, since they 
represent the departure from an exact initial data set (see □ for their definitions) . 
In this example, we excise the points inside a sphere of radius a + h centered at 
the hole, where the absolute value of the specific angular momentum a is also the 
radius of the Kerr ring-like singularity and h is the grid spacing (for our example, 
h = M/4). Figure [j] shows the Hamiltonian constraint residuals near the rightmost 
hole of the head-on collision. The circles show the second order truncation error 
and the squares the residuals from the attenuation method. First we note that in 
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Figure 1. Hamiltonian constraint residuals near the rightmost hole of the head-on collision. 

some areas the violation of the constraints for the attenuated data is greater than 
that corresponding to truncation error. However, even in these areas the viola- 
tion is smooth and small in absolute value, as opposed to the unbounded behavior 
shown by the truncation error near the singularity: the divergent behavior present 
at the singularity has disappeared with the use of attenuated data. This is due 
to the fact that, having cancelled the influence of the presence of the second hole, 
the metric and extrinsic curvature become the fields corresponding to an isolated 
black hole at the location of the singularity. This simplifies the exact problem, since 
this method provides exact inner boundary conditions for the elliptic solver. The 
applicability of these concepts to the numerical solution of the initial value problem 
has been confirmed by Marronetti and Matzner El, who presented the first results 
for a black-hole binary using the attenuation method to generate background data. 

These methods can be further refined to provide more astrophysically realistic 
initial data for the case of black holes in circular orbits, by deriving the extrinsic 
curvature Kij from the presence of a Killing vector £ = dt + t-ud^, instead of using 
the superposition of two boosted holes. Work under progress will incorporate these 
variations into the full numerical problem. 
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